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the cost of moving extra power in the 
network, the penalties incurred for 
late train departures due to the failure 
to assign an appropriately-equipped 
locomotive and the cost of train 
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wo 01/49548 PCTAJSOl/00119 



OPTIMAL LOCOMOTIVE ASSIGNMENT FOR A RAILROAD NETWORK 

5 FIELD OF THE INVENTION 

This invention relates to a method for determining the optimal assignment 
(distribution) of locomotives on a raibroad network, by taking into account various 
costs and penalty factors associated with the movement and distribution of the 
locomotives. 

10 

BACKGROUND OF THE INVENTION 

One problem faced by raibroad operators is the assignment of a sufficient 
number of locomotives of different types to operate on a pre-planned schedule. For 
efficient and timely operation of a railroad, it is essential tiiat the tractive power (i.e., 

15 locomotives) is distributed around the network as required to move people and 
freight. The required distribution is also dynamic with respect to time and constantiy 
changing due to track outages, planned or emergency locomotive maintenance, 
weather conditions, etc. 

The problem of locomotive assignment can be illustrated by the following 

20 simple example. A railroad moves grain from Nebraska to a port on the Gulf of 
Mexico. To move the grain south requires a certain amount of tractive effort provided 
by the railroad locomotives. However, moving the empty cars north back to 
Nebraska, requires far less locomotive power because the cars are empty. The result 
is an accimiulation of tractive power at the Gulf of Mexico port. To minimize its 

25 costs, the raibroad desires to move the locomotives back north to Nebraska in a cost 
efficient maimer. For example, the locomotives could be put into service moving 
freight on a different route such that they eventually find flieir way back to Nebraska, 
where they can again be put into service hauUng grain south. 

Each railroad has at least one locomotive dispatcher with the responsibility of 

30 ensuring that each terminal has the correct number of locomotives to move the freight 
on schedule. The locomotives must be in the yard and ready to depart on time. To 
carry out this fimction, the dispatcher must assess the location, movement, and 
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availability of each locomotive and then predict when and where it will be available 
for the next train consist. For example, the dispatcher may be operating 16 hours 
ahead in planning the movement and availability of locomotives. 

The cost of moving a locomotive either with a train as active power, with a 
5 train in-tow, or by itself must be considered by the dispatcher. A penalty may be 
incurred by the railroad, under its contract with a customer if a train is not ready to 
depart on time because there is not sufficient tractive power to move it or if the train 
arrives late at its destmation. The penalty applied to a late arriving train is actually a 
sum of penalties associated with each individual car. The penalty incurred for each 
10 car depends on the shipping commitment pertinent to that car. The penalty could be 
zero for cars that have arrived ahead of schedule and could be very large for cars that 
are significantly behind schedule. In some situations, it may actually be cheaper for 
the railroad to pay the late penalty than to incur the cost of moving a locomotive to 
meet the schedule. 

15 The dispatcher must also give consideration to various circumstances that may 

interrupt the smooth flow of freight on the railroad network. The dispatcher cannot 
exceed the locomotive capacity of each yard, and the traffic on each railroad corridor 
must be established to ensure that at any given time there are not an excess number of 
trains travelling the corridor, which would slow down the operation of all trains on the 

20 corridor. Locomotives must undergo various levels of inspections at predetermined 
intervals. Sunple mspections are performed in the yard, while more thorough 
inspections must be performed in an inspection shop. The timing and duration of 
these inspections must be accounted for in determining the optimum distribution of 
locomotives in the network. Further, the productivity of each inspection shop will 

25 determine the length of time the locomotive is out of service. 

Locomotives are segregated into power classes, with each power class having 
a different power/speed relationship, that is for a given speed the locomotive can 
generate a given amount of tractive power. The tractive force capacity influences the 
train consists to which the locomotive can be assigned. Obviously, heavier trains 

30 require the use of locomotives capable of pulling heavier loads. Assignment of 
locomotives must also be accomplished with consideration to track rail and load 
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bearing capacities. Locomotives that are too heavy can damage the rail or the 
roadbed. 

Raihoad dispatchers communicate with the locomotive engineers by sending 
signals through the rails and through over-the-air communications links. There are 
two popular signahng systems referred to as Centralized Train Control (CTC) and 
Automated Block Signaling (ABS). For both the CTC and ABS signaling systems the 
track is divided into sections or segments, with each segment having a detection 
circuit. When a train enters a section, it closes the electrical circmt and sensors 
identify that a train occupies the section. Centralized Train Control comprises a 
system of computers that read the track sensor status and transmit data regarding train 
positions to a central dispatcher. The central dispatcher (either a computer or a 
human operator using a computer) schedules and controls train movements by setting 
the color of each individual traffic signal based on train position and specific railroad 
operational rules. In a centralized train control system, all signals in an area are 
controlled firom one central point In the ABS system, traffic signals colors are 
controlled by the train on the track. The color of the traffic signal at the beginning of 
the occupied segment is red, the color of the traffic signal on the previous segment is 
yellow, and the color of the signal on the next previous segment is green. Thus, the 
color of the signal is determined from the status of a single track segment. 

In places along the railroad track where there is reduced visibility (narrow 
curves, steep slopes, etc.) secondary signals are installed. The secondary signals 
"repeat" or "predict" the color of the next primary signal using combinations of colors 
and light positions. In the ABS system, a secondary signal is always subservient to its 
main signal, while in the CTC system each signal is independently controlled and, 
therefore, can change its aspect. Both the ABS and CTC systems send signals 
identifying the status of the next primary signal through the rails to the locomotive. 
Special equipment installed in the locomotive cab receives the signals to advise the 
engineer of the color of the next traffic signal before it becomes visible. 
Unfortunately, the systems use a different code of electrical signals for transmitting 
the signaling information. 

Some locomotives are equipped to respond only to CTC signaling, others only 
to ABS signaling, and some locomotives to both. Occasionally, a locomotive may 
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have no on-board signaling capability. Because the signaling system advises the 
engineer as to the conditions of the rail, occupancy of subsequent rail blocks along the 
path of traivel, etc., the lead locomotive must have a signaling system capable of 
interfacing with the signaling system of the track over which it travels. 
5 As can be seen, the schediiling and dispatching of locomotives is an extremely 

complex problem that can cause significant inefficiencies and extra costs to the 
railroad. In many situations, there simply is no solution to an under capacity problem 
other than rmming the locomotive; by itself, to the place where it is needed. This is an 
obvious cost to the railroad in terms of personnel, fuel, and wear and tear on the 
10 locomotive. Then there is the additional problem of finding a time slot in the corridor 
schedule so that the locomotive can travel to the temiinal where it will later be 
needed. 

SUMMARY OF THE INVENTION 

15 In accordance with the present invention, the above described shortcomings of 

the locomotive dispatching process are obviated by a new and improved dispatching 
algorithm. The process utilizes a mathematical model of tiie locomotive assignment 
problem that assigns the locomotives in the network to various terminals at minimum 
cost and at the appropriate time. The present invention minimizes the costs of moving 

20 locomotives to achieve the desired distribution at the appropriate time and minimizes 
the penalties for delaying train departures. The present invention creates a plan for 
the distribution of all locomotives owned or controlled by a railroad with regard to 
power class and the cab signaling capabilities of each locomotive. The process can be 
executed using a sliding time window of any duration fi:om the current moment to the 

25 future. 

The present invention plans the movement of locomotives across the railroad 
network as either an active locomotive or as extra power connected to another train 
consist (i.e., in tow), and accomplishes optimal dispatching by considering the various 
parameters influencing the use and movement of locomotives across the network. 
30 These parameters include the power class and signaling capabilities of each 
locomotive. The process takes into account the number of locomotives needed at 
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10 



15 



20 



25 



various points in the network and various costs associated with moving or holding the 
locomotives so that they are available when needed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention can be more easily understood, and the further 
advantages and uses thereof more readily apparent, when considered in view of the 
description of the preferred embodiments and the following figures in which: 

Figure 1 is a block diagram of a locomotive movement planner including the 
schedule optimizer of the present invention; 

Figure 2 is a flow diagram illustrating operation of the present invention; and 

Figure 3 is a block diagram illustrating the incorporation of the present 
invention into a real-time railroad dispatching system. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention is a process for optimizing the assigmnent of 
locomotives in a raiboad network. The process performs the planning of locomotive 
assigranents across a selected area, and minimizes the cost of locomotive distribution, 
the penalties paid for delaying train departures and the cost of performing the 
operation of moving locomotives to the point where they are needed at the appropriate 
time. The present invention creates a plan for the distribution of aU locomotives 
owned, leased or controlled by the railroad, with regard to the power class and cab- 
signaling capabilities of each locomotive. The process can be executed using a 
sliding time window of any size, firom the current moment to several days (or longer) 
into the future. The process plans the movement of locomotives across the network as 
active power (providing tractive power to a consist) or as extra power (in tow coupled 
to another locomotive). 

Definitions and Notations 

All known mathematical values are printed in upper case, while unknown 
variables are in lower case. Sets are printed in upper case bold character (for 
example, !). 
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Each locomotive operating on the railroad network has two descriptors. The 
first descriptor represents the locomotive power class (C36-7, SD45-T2, C44AC, 
etc). Since the power class is a part of the technical specifications of the locomotive, 
it cannot change. In the present invention locomotives are defined by their power 
5 group, which accounts for the limited compatibility that exists between certain power 
classes. Locomotives in the same power group can be coupled together. The process 
includes a matrix defining ttie compatible classes of power. 

The second descriptor defines the cab signaling system installed on the 
locomotive. There are two cab signaling systems in use on North American railroads 
10 today. These are the Centralized Train Control and the Absolute Block Signaling. As 
discussed above, some locomotives are equipped with only the CTC system and 
others are equipped with only the ABS system. Other locomotives may not have a 
signaling system capability, while still oth^s may have both the ABS and CTC 
systems installed. 

15 The cab-signaling category changes if the cab equipment is replaced. Note that 

the cab-signaling category is important only for the lead locomotive (typically, a first 
locomotive at the firont of the train), firom where the engmeer drives the train. Any 
locomotive that cannot be utilized as a lead locomotive, for example because it is not 
winterized during the cold season, has cracked windows, or has faulty non-vital 

20 control equipment, is assigned a signaling category "none" to account for its inability 
to serve as the lead locomotive. 

The square matrix S's provides the mapping between the signaling 
equipment installed in the locomotive cab and the signaling equipment installed on 
the corridor over which the locomotive is traveling. The upper index refers to flxe 

25 signalling system installed on the track, and the lower index refers to the cab-installed 
signalling equipment. The values on the main diagonal of the matrix are equal to one, 
indicating that the cab signaling equipment is identical to the track installed 
equipment. The exception is the line and the column associated with "none" or no 
signaling equipment installed, where the entries are zero. 

30 The matrix values for track-installed signaling systems that are a subset of the 

locomotive-installed equipment are also set to one. For example, the compatibility 
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between a locomotive having both CTC and ABS equipment and the track signaling 

CTC 

system CTC is Sctc+abs"^^- The value of one shows the fact that locomotives 
equipped with CTC and ABS signaling can be used as lead locomotive for trips 

requiring CTC signaling only. Needless to say, Octc+abs^^ too. Obviously, 

5 OcTc ~^ since locomotives with only CTC installed system cannot cover track 
corridors with CTC on one section and ABS on the other section. 

The present invention includes some needed redundancy. Some parameters 
refer to a train t originating at location i. The redundancy comes from the fact that a 
train can originate at only one location, and, once the train is identified, the 
10 originating location becomes unambiguous. But, two indices are used to describe the 
train because some operations in the process require summation over either one of the 
two indices. 

The term "power change location" as used herein, designates a track length 
that is small enough to ignore the time and cost for moving locomotives along it. 

15 Such a section of track could be as small as a yard or a terminal, but can also be much 
larger, as will be shown below. The process performs the planning for moving power 
between such locations, but establishing the extent of each track length is a matter of 
implementation and may depend of many factors that are specific to each situation. 
Keeping the niunber of such "power change locations*' as small as possible reduces 

20 the amount of data on which the process must operate and thus will decrease its run 
time. 

Table 1 below presents six categories of entities: indices and parameters, sets, 
unknown, static, dynamic and auxiliary entities, which are used in the process of the 
present invention. Indices and parameters are used for qualifying entities used in the 

25 system of constraints. All entities quaUfied by indices and parameters are organized 
in sets: sets of locomotive, sets of trains moving between two locations, etc. The 
static entities are not time dependent and are not modified during the optimization 
process performed by the present invention. The value or values of the static entities 
are well known before the optimization starts. For example, the amount of power 

30 generated by a locomotive is independent of the optimization process. The penalty 
paid per time unit for not powering a train is known before the optimization process 
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starts, although it may have different values from one time period to another. For 
example, it may be zero for the first few hours, when all of the cars will be able to 
make their connection to next train; after that the penalty may jump as a block of cars 
DMSS the next connection, jeopardizing on time deUvery. 
5 Values of dynamic entities are known at the beginning of the optimization 

process, but for the reminder of the time window of interest they depend on values of 
certain unknown variables that are associated with previous stages of the process. For 
example, the number of locomotives at one particular location is a known entity when 
the optimization starts. The optimization process programs the movement of 
10 locomotives across the network, which makes the entity, "number of locomotives at 
the given location" dependent of elements whose values are found while optimizing 
algorithm runs. 

Auxiliary entities are used only for making the presentation easier to 
imderstand. 

15 Table 1 



Symbol 


Meaning 


Comments 


Indices and parameters 


C 


Indices for differentiating the class of power. 




G 


Indices for marking the locomotive group of 
power. 




hj 


Indices indicating one location or two adjacent 
locations (with direct connection). 




S 


Indices indicating the locomotive cab signaling 
category. 




r 


The time parameter. 


Integer values. 


Sets 


C 


The set of classes of power. 




G 


The set of groups of power. 




K 


The set of corridors in network. 




Nj 


The set of locations adjacent to location 




S 


The set of signaling categories. 
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Symbol 


Meaning 


Comments 




The set of trains at location z planned to move 
to location / 




Y 


The set of locations in network. 




Unknown entities 


a,: (^) 


J. ill? ilUllil^d \jx IvUvliXUU Vwo jJ\JVv\^x. y^xaoa C' \^au 

signaling type s assigned to train t initiated at 
location i during the time interval {v, v^l). 


Jxiic^cJi vcuuwd xjuxy. 




Indicator for marking that train / at location i is 
powered with locomotives from power group g. 


Integer values only. From the 
system of constraints will 
result that it can be 1 for no 
more than one power group. 




Number of locomotives of signaling category 5* 
power class c that are on hand to be assigned 
to trains in location i during the time interval 


Integer values only. 




Number of locomotives of signaling catecorViS 
power class c that arrive as active power at 
location i during the time interval (r, v^l). 






Confirms the utilization of crew planned for 
train t originating at location i during liie time 

interval (r,r4-7/ 


Values of zero or one only. 




Number of locomotives of signaling category 
5, power class c, planned to exit location i and 
directed to adjacent location j as extra power 
during the time interval (r, v^l). 


This entity does not include 
the active power that actually 
pulls the train. The extra 
power is moved from one 
location to another in tow. 




Number of locomotives of signaling category 
s, power class c arriving at location i from 
adjacent location J as extra power during the 
time interval (z; r+7j. 


Same as above. 
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Symbol 


Meaning 


Comments 


Static entities 




Number of locomotive of class c and signaling 
category s planned to arrive at location i from 
direction j during the time interval (r, Tfi). 


These are locomotives 
already on route when the 
optimization begins. 


C 


The cost of using a locomotive of class c as 

active power at maximum capacity. 


In $ per mile. 


Ki.j 


The capacity of the corridor between 
locations i and J. 




A., 


The traveling distance of train t originating at 
location /. 


In miles. 




The distance between location i and adjacent 
location j\ 


In miles. 


r 


The cost of moving a locomotive of class c as 
extra power. 


In $ per mile. 




The matrix shows the equivalence 
(compatibility) of different classes of power. 


Values 0 (not compatible) 
or 1 (compatible). 




The upper limit of yard capacity. 


Ninnber of trains. 


Mi,, 


Maximum number of locomotives that can be 
attached to train t originating at location i. 


Prevents assigning of too 
many small locomotives to 

a very heavy train. 


a. T 


Minimum number of locomotives that can be 
attached to train t originating at location i. 


Prevents assigning too few 
locomotives to a train. 


n;(*;Ar-f 


The probabiUty that a locomotive of power 
class c will be released from maintenance at 
location i between k and k+1 hours. 


Values in this table reflect 
technical capabilities of the 
repair shop at location i as 
well the frequency of 
maintenance and repair 
operations. 


) 


Power required at location / for train t using 
signaling category so during time interval 


In 1000 HP. The model 
considers that the power is 
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Symbol 


Meaning 


Comments 




(r,zH-i). 


available and can be 
scheduled to move to the 
departure yard within the 
{r,T\rl) interval. 


s': 


Signaling compatibility matrix. 


Equal to one if signaling 
system s installed on 
locomotive can support the 
track signaling requirement 
so. 


<S>,j 


The time for traveling from location / to 
adjacent location j. 


In time intervals. 




Time needed to complete the departure 
procedures at location i (travel from locomotive 
shop to departure yard, charging the brake 
system, inspection, brake testing, etc.). 


In time intervals. 


ir 


The power of one locomotive in power class 
c. 


In 1000 HP. 


Dynamic entities 




Number of locomotives of power class c and 
cab signaling s available at location i at the 

beginning of the time interval (r, v^l). 


Integer number. 


AW 


Nimiber of trains at location i at the begirming 
of the time interval (r, v^l). 






Number of locomotives of power class c and 
cab signaling s entering the fleet at location i 
during the time interval (t, t^I). 


These locomotives can be 
on loan, new acquisitions, 
refurbished or foreign. 




Penalty for not providing power to train t at 
location i during time interval (r, v+l). 


In $ per time period. 

/>^^(T;r < To) = 0 where tq 

is the time when the train 
power should be scheduled. 
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Symbol 


Meaning 


Comments 




Number of locomotives of power class c and 
cab signaling s exiting the fleet at location i 
during the time interval (r, z+l). 


Those locomotives can 
retire, be removed from the 
fleet for major repairs or be 
returned to a foreign owner. 


Yr la V.*'' 


The availability plan for the crew operating the 
lead locomotive of class c on train / originating 
at location i during the time interval (r, v¥l). 


Plan generated by the crew 
management system. 




Number of locomotives of power class c and cab 
signaling s released from maintenance shops at 
location i during the time interval (t, t^I). 




Auxiliary entities 




Logical variable associated with train t at 
location i during the time interval (r, T+i/ 


1 = train gets reqidred 
power during the time 
interval (j, t^1)» 
0 = train does not get power 
during the time interval 




Number of locomotives of category s and 
power class c assigned to trains leaving the 
location / during the time interval (r, r+l). 





The Algorithmic Process 

The process includes seven elements: local power assigEunent, network 
balance, receiving yard balance, locomotive shop release, network-location balance, 
5 dynamic relations and the objective ftinction. 

The local power assignment element deals with the operation, at the location 
of interest, of assigning locomotives to trains. The network balance contains 
constraints linking each location with other locations within the network. The 
receiving yard balance counts the number of locomotives arriving at the location 
10 within each time interval. The locomotive shop operations are represented as delays 
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between the moment a locomotive enters the receiving yard until it becomes available 
for assignment. The network-location constraints assure that the number of 
locomotives assigned to all locations is not larger than the number of locomotives 
actually available in the system. In another embodiment, this section could reflect the 
rental availability of locomotives. The dynamic relations link one time interval with 
the next time interval. 

All these elements are represented by a system of constraints that can have 
many feasible solutions. The best solution is the one that minimizes: the cost of 
moving extra power across the network, the total penalties paid on the whole network 
and the total direct cost of moving the freight across the network. The process 
presented is in a class of mixed integer linear problems that can be solved using either 
linear methods (known as Danzig- Wolfe methods) or interior point methods. Such 
minimization solution methods are well known by those skilled in the art. 

Local power assignment. 

The assignment of locomotive to trains should reflect the following facts: 

a) a train should either have as much power as required, or none (that is, 
the train is delayed or cancelled), 

b) the set of locomotives assigned to a tram should be from the same 
power group, 

c) the set of locomotives assigned to train should have at least one 
locomotive (used as the lead) supporting the required cab signaling 
system of the track over which it will travel. 

seS ceC 

alMh Z* VieTj e U T,j,ge O (3-1) 

< To -» 0 

Equation (3-1) assures that whatever locomotive assignment is made, the 
assigned power is at least equal to the required power. According to this equation, the 
at W variables can have only integer, not negative, values. 
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Eaf,(^)=^M(^)-^».(^-i) VieF.r e U T,j (3-2) 
/ ^ . ^ V/ey,/€ U T,, (3-3) 

Equation 3-2 accounts for the dynamic characteristics. According to equation 
(3-2), all values of the a?,W variables associated with the same train and.representing 
5 all power groups, added together, equal no more than one. Since all a'M) are non- 
negative integers, this means that no more than one value of the variables is equal to 
one, wliile all others are zero. This reflects the fact that all assigned locomotives to the 
same train should be from the same power group. 

The train travels in a region requiring signaling capability of category so- This 
10 means that from all assigned locomotives at least one should support this category of 
signaling. This locomotive will occupy the lead position in the train. 

I^E'-'HS^a^i^^aU^) v.-^ir,re U T,,.^eG (3-4) 

ceC sea ^^"f 

The inequality sign used in equation (3-4) makes sure that at least one 
locomotive of the assigned power group has the equipment for signaling category 5, 
15 but assigning more than one such locomotive is acceptable. 

The crew driving the train should be qualified to operate a locomotive with 
power class c and with the itinerary of train / originating at location i. 

Equation (3-5) assures that the locomotive selected as lead can be operated by 
20 the crew mentioned in the availability plan. Since federal regulations require that 
engineers be qualified for the trip route, not necessarily for the type of equipment 
used on the trip, all values of W]^l^^) should be equal across the c coordinate in the 

United States. That is, the values within the matrix will be identical along the 
dimension c (for each value of c). There will be no differentiation caused by different 
25 values of index c specifying the type of equipment. In other countries where there are 
many different types of locomotives, each requiring a particular skill set to operate, or 
where the engineers are required to perform maintenance on the engine, the situation 
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would be expressed differently by the matrix. A particular engineer might be 
qualified for only some type of equipment, which would change the availability 
matrix along the dimension c. In other words, for different values of c the matrix 
could have different values, even if all other indices are kept constant. 
5 The total number of locomotives used for powering all trains at location / by 

each power class c and cab-signaling s is computed with equation (3-6). 

E a;;(^)=wr(0 V,er.5eS.ceC (3-6) 

Equation (3-7) restricts the number of locomotives that can be used for 
powering a train. 

10 N,,^YLaTM)^Mi. V,-eF./eUT,, (3-7) 

The relation has the effect of preventing the assignment of too many or too 
few locomotives to a heavy train. 

Network balance 

15 Equations in this category are responsible for balancing the locomotive power 

in the network. The number of locomotives for cab signaling category s and power 

class c available at location i during the time interval (t, t^I) is p^"" (r) 

yi^Y,s^S,c^C (3-8) 

20 Equation (3-8) indicates the total number of locomotives available to be 

assigned at location i, taking into account the locomotives already available at the 
location, the inconaing extra power from all neighboring locations, the extra power 
leaving for all neighboring locations, locomotives released from maintenance shops, 
as well as locomotives entering and exiting the fleet at location i . 

25 The extra power arriving in a yard during the time interval (t,t^1) left the 

originating location **hours" ago, as presented in equation (3-9). The extra 

power should be coming in tow to location /, therefore it does will not need any 
maintenance and can be immediately assigned to trains departing from location / . 

15 
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J^^i Vi,7 G e S,c e C (3-9) 

x'j{r;T <0)=0 



Receiving Yard Balance 

The number of locomotives arriving at location / as active power during the 
5 time interval (%r+l) depends on the departure time from other locations and the 
duration of each trip. 

Pr(0=ZZ«5:?(^-6>y.-(9>i:^;;?(^) V.eF,,eS.ceC (3-10) 

Equation (3-10) shows that locomotives arriving within the time interval 

(t, T^l) have been subject to an assignment at time "^"Q^^Ou • 

10 Locomotive shop release 

Locomotives are released from the maintenance shop according to the length 
of the maintenance procedure and/or inspection and shop productivity. 



L * 



V/ey,jeS,ceC (3-11) 



The expression between brackets should be rounded to the nearest integer to 
15 provide the number of locomotives. 

Connection Between Network Availability and Local Request 

Obviously, the total number of locomotives used at a location defined in 
equation (3-6), caimot be larger than the nimiber available at the location as defined in 
20 equation (3-8). This is expressed below in equation (3-12). 

uT{r)< p:\r) \fi e T,s e S,c e C (3-12) 
Corridor Capacity 

The nimiber of trains scheduled over a corridor should not be larger than the 
25 corridor capacity for time interval of interest. 
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'ZlKi^)^^>j "^i^rjsN, (3-13) 

This relation restricts the number of trains leaving from location i toward 
location j to the specified number of trains per planning interval. 

5 Yard Capacity 

As expressed below, the total number of trains in a yard at any particular time 
should be smaller than yard capacity. 



^^(^■''^ V/^y (3-14) 



10 Dynamic Relations 

Now we consider the set of equations making the link between variables 
representing one time interval and variables of the next time interval. 

X(r + 1) = p;\T)-u;{r) Vier.56S,ceC (3-15) 

Equation (3-15) computes the availabiUty of locomotives at the be g i nnin g of a 
15 new time interval as the difference between the total power available during the 
previous time interval, minus those locomotives that were assigned to trains. Based on 

equation (3-12), A't'^i^ + 0 is always a non-negative number. 

Because the time when a train can depart is known in advance, the S3rstem of 
equations should be written so that the train can depart at anytime, and then let the 
20 optimization process discussed below select the best departure time. 

Equation (3-2) included the variable ^i,ti^). Its value is one only if the train 
has enough power and can leave the yard. 

Both equations labeled (3-16) above (note, the second equation of (3-16) is 
25 simply a rewrite of equation (3-3) at time (z+7) propagate the value of the ^lA^) 
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variable from the moment the train was assigned power, forward in time. An example 
of the time sequence of values for this variable for the same train is: 
OA'-'jOjljl*!*-- *!* where the first zero corresponds to the planned time when the 
power should be assign (that is, the estimated time of departure less the duration of 
5 departure yard operations), and the first one corresponds to the time when the train 
was assigned power. 

Because S^X^) acts like a step function, equation (3-2) provides a set of 

constraints where only one value of ^^^^ (t) can be equal to one at any time. It assures 
that all locomotives assigned to a train are in one power group only. 

10 

Objective Function 

The objective fimction below provides a measure of profitability of the whole 
operation. To perform a truly economical optimization, it would be necessary to 
know the unit cost and unit revenue for each type of service provided, plus general 

15 overhead costs. With this additional information, a complete economic analysis could 
be performed and the business profitability optimized. The optimization analysis here 
is therefore incomplete because there is no economic analysis. The optimization 
process does not necessarily imply that the operation will be more profitable if the 
derived optimization solution is implemented. The implementation will no doubt cut 

20 costs, but may, in fact, also create a revenue decline. 

Generally, with no additional constraints and using a criteria of cost 
minimization only, in the limit the conclusion would be to close a business because 
that is the only way to achieve the smallest possible cost (i.e., zero). The algorithm of 
the present invention attempts to minimize the cost using a constraint that all goods 

25 are to be moved in accord with the customer's schedxile. In this algorithm we 
minimize the global cost of the operation. 
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min 



fsF JeN/seQ cgC 



/er te Jr. 



l,J 



z z zza;;wc'A 

/ey <e IJr, y *€S ceC 



(3-17) 



The first term under the summation evaluates the cost of movmg the extra 
power in the network. It is the product of the nximber of moved locomotives, the 
distance between locations and the cost per mile for towing a locomotive. 
5 The second tenn evaluates the penalties paid for not powering a train in time. 

The 5^ f (r) variable has the value one when the train has power to move and zero 

when it cannot move. The objective ftinction has a summation over time, which 
makes it possible to count the cost of departure delays. The value of the penalty per 
time unit is static, but it can be different from one hour to another. 
10 The third term evaluates the cost of operating the train using a particular 

power class. This term tells the optimization process to select the most productive 

equipment for doing the job. The cost is in proportion to the distance Ai^ and the cost 

per mile C*^ associated with the power class. 

Those skilled in the optimization art will realize that the size of the model can 
15 be reduced using algebraic substitutions. Such reduction is useful for decreasing the 
size of the data involved in the optimization process, therefore increasing the 
execution speed. 

Equation (3-1) is used in the form 

ZZa:;(^)C/'£''"-af.Wi2::(0^o v/er./^ yT,,.^eG (5-1) 

20 Equation (5-1) will generate |T||G( constraints, where |X| is the cardinality of 

setX. 

Equations (3-2), (3-3) and (3-4) are modified as follows, in one embodiment 
of the present invention. 



19 



wo 01/49548 



PCT/USOl/00119 



ZalM-^M^ ^.A^ - 1)= 0 Vi e T,t e U T,, (5-2) 
S„{t)<1 \fieY,tG U T,j (5-3) 

i:E:"i:s':d;;i')-ai{->^ yi^y,t^}j%,,g^o (5.3) 

The number of constraints generated from equations (5-2) and (5-3) is equal to 
5 2|T|, while equation (5-3) generates |T1|G| constraints. If crew management is not of 
interest, equation (3-5) should not be used. 

The Tnini-murn ' limit on number of locomotives attached to a train is not a 
concern. In such a case, equation (3-7) can be used for limiting the maximum number 
of trains only. Equation (5-5) below, derived from equation (3-7), generates |T1 
10 constraints. 

ZZ^ZnC^)^^.-. V/er,feL)T,, (5-5) 

The network balance equation (3-8) can be combined with (3-6), (3-9), (3-10), 
(3-1 1) and (3-15). The result is: 



15 Vi€F,5€S,ceC (5-6) 

The network balance is considered for all yards, for each signaling category 
and for each power class. Equation (5-6) generates |Y||S||C| constraints. 
The constraints for corridor and yard capacity are 



20 



wo 01/49548 



PCT/USOl/00119 



10 



15 



The total number of constraints generated witii equation (5-7) is |K|, the 

number of corridors in network. 

r \ 



1^1 



-/l(r + l)=0 



VieF (5-8) 



Equations (5-8) generate 2|Y| more constraints. 

The train assignment step function is controlled by the relation 
<y,,,(T +!)-<?,,, (t)^0 VfeF,/6UT-.7 (5-9) 

It generates another |T| constraints. 

The second term of the objective function is decomposed and the constant 
value 5] 2 2 Pih^ ^™ ignored while searching for the optimal solution. 



The function to optimize is then: 

( r 



mm 



s z zz^;;wcrA 



(5-10) 



The total number of constraints for this problem is then 

2\T\\G\^4\Ti+\¥\mC\+\K\^2\¥\ (5-1 1) 

The nimiber of variables is computed as follows: 

a;;(r) |T||S||C| 

a'M |T||G| 

|T| 

Xx?;(t) |k||s||C| 
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The linear problem generated with this model has a number of slack variables 
(compensating for the inequaUties) equal to 2|T||G|+3|T|+|K|+|Y|. In order to find the 
size of the problem, the number of imknown variables, slack variables and constraints 
are multipUed by the number of time intervals considered in problem horizon. Those 
5 skilled in the art will reaUze that various computer-based techniques can be employed 
for optimizing the objective function, including the techniques set forth in Large Scale 
Linear and Integer Opt imization- A Unified Approach, by Richard Kipp Martin, 1999. 

Figure 1 is a block diagram of a locomotive assignment planner 10 constructed 
according to the teachings of the present invention. The locomotive assigmnent 

10 planner 10 includes an optimizer 12, incorporating the process of the present 
invention, that receives as an input the stored train schedule firom an element 14 and 
the current status of all locomotives firom an element 16. The stored train schedule 
includes the desired departure and arrival time of the trains operating on the network. 
The output fi-om the optimizer 12 is the optinaized locomotive assigmnents, 

15 The flow chart of Figure 2 illustrates the operation of the present invention. 

At a step 22, the constraints set forth in equations (5-1) tiurough (5-9) are identified. 
The objective fimction of equation (5-10) is calculated at a step 24. At a step 26 the 
objective function is minimized to provide the optimized locomotive assigmnents. As 
is well known by those skilled in the art, these steps can be carried out by a general 

20 purpose computer. 

Figure 3 is a block diagram of the real-time dispatcher system 30. The 
locomotive assignment planner 10 provides the optimized locomotive assigmnent 
schedule to a real-time control system 32, which dispatches the trains and controls the 
track switches. The real-time control system 32 is in communications with 

25 locomotives 34 via a communications device 36. Dispatching information can be sent 
to the locomotives 34 either through an RF Unk or through a track-based system. The 
real-time control system 32 also controls the track switches (illustrated by element 38) 
to properly route the trains, via the conmiunications device 36, 
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WHAT IS CLAIMED IS: 

1. A method for optimizing the assignment of locomotives in a raikoad 
network by roinimizing the cost of moving locomotives in the network over a 
predetermined time period, r , wherein each locomotive is within a signaling class s 

5 from a set of signaling classes S, and wherein each locomotive is within a power class 
c from a set of power classes C, and wherein the distance D^ j between locations i and 

J on the network is known, and the cost of moving locomotives as extra power as 
a function of power class is known, said method comprising the steps of: 

(a) determining the number of locomotives of signaling class s and power 
10 class c that are planned to exit location i in the direction of location j as extra power 

during the predetermined time period r ; 

(b) determining the product of the result of step (a) for each power class, 
the distance between locations i and j\ and the cost of moving locomotives as extra 
power as a function of each power class c; 

15 (c) summing the result of step (b) over all power classes; 

(d) summing the result of step (c) over all signaling classes; 

(e) simiming the result of step (d) over all adjacent locationsy; 

(f) summing the result of step (e) over all locations i; 

(g) sunaming the result of step (f) over the predetermined time period r ; 

20 and 

(h) minimizing the result of step (g). 

2. The method of claim 1 wherein the execution of step (h) is subject to 
the following constraint, 

sea ceC ■'^ ' 

25 3. The method of claim 1 wherein the execution of step (h) is subject to 

the following constraints, 

ZaM'^)- '^/.(^ -1)= 0 Vi e T,t e T,, 

ceC leS J^^i 
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4. The method of claim 1 wherein the execution of step (h) is subject to 
the following constraint, 

5. The method of claim 1 wherein the execution of step (h) is subject to 
the following constraints. 



/J 



nf(k)>of 



nf{k) 



6. The method of claim 1 wherein the execution of step (h) is subject to 
the following constraint, 

7. The method of claim 1 wherein the execution of step (h) is subject to 
the following constraints, 



geG 



f 



-A(r + l)=0 



A,(r)<Z,. 

8. The method of claim 1 wherein the execution of step (h) is subject to 
the following constraint, 

9. The method of claim 1 including additional step (i) of scheduling the 
movement of the locomotives based on the result of step (h). 
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10. A method for optimizing the assignment of locomotives in a raikoad 
network by minhiiizing the penalties associated with late train departures, wherein the 
penalty incurred Pij{t) if a train f at a location i does not have power over the 

predetermined time period r is known, said method comprising the steps of: 
5 (a) determining whether a train t is powered during the predetermined time 

interval r ; 

(b) assigning a first value to each train t if the train t does not have power 
and assigning a second value to each train t if the train / does have power; 

(c) multiplying tiie result of step (b) by the penalty incurred; 

10 (d) summing the result of step (c) over all trains t moving between 

locations i and j; 

(e) summing the result of step (d) over all locations i; 

(f) summing the result of step (e) over the predetermined time period r ; 

and 

15 (g) minimizing the result of step (f). 

11. The method of claim 10 wherein the execution of step (g) is subject to 
the following constraint, 

12. The method of claim 10 wherein the execution of step (g) is subject to 

20 the following constraints, 

d,,{r)<\ Vier,/e U T,, 

Sr-^E^r d;:i^)-al{r)^^ e Y,t e U %„g e G . 

13. The method of claim 10 wherein the execution of step (g) is subject to 
25 the following constraint, 

14. The method of claim 10 wherein the execution of step (g) is subject to 
the following constraints, 

25 



wo 01/49548 PCT/USOl/00119 



5 



10 



15 



nf{k)>o 



k=0 



15. The method of claim 10 wherein the execution of step (g) is subject to 
the following constraint, 

16, The method of claim 10 wherein the execution of step (g) is subject to 
the following constraints, 

r \ 

-yl(r + l)=0 



geO 



4(r)<i, 

17. The method of claim 10 wherein the execution of step (g) is subject to 
the following constraint, 

ye-"! 

18. The method of claim 10 including the additional step (h) scheduling 
the movement of the locomotives based on the result of step (g). 

19. A method for optimizing the assignment of locomotives in a railroad 
network by minimizing the cost of operating trains using a power class c, wherein the 
cost of using a locomotive of power class c as active power is known and wherein the 
traveling distance of each train departing from location i is known, said method 
comprising the steps of: 
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(a) determming the number of locomotives of power class c and signaling 
class s assigned to a train t initiated at location i over the predetermined time r ; 

(b) calculating the product of the result of step (a), the cost of using a 
locomotive of power class c as active power in a train t and the distance train t will 

5 travel; 

(c) summing the result of step (b) over all power classes; 

(d) summing the result of step (c) over all signaling classes; 

(e) summing the result of step (d) over all trains traveling between 
locations z and j; 

10 (f) summing the result of step (e) over all locations i on the network; 

(g) summing the result of step (f) over the predetermined time r ; and 

(h) minimizing flie result of step (g). 

20. The method of claim 19 wherein the execution of step (h) is subject to 
the following constraint, 

5€S CeC ' 

21. The method of claim 19 wherein the execution of step (h) is subject to 
the following constraints, 

<J),,(r)^l V/ey.^e U T,, 

22. The method of claim 19 wherein the execution of step (h) is subject to 
the following constraint, 

23. The method of claim 19 wherein the execution of step (h) is subject to 
25 the following constraints. 
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JeN, JeNf 



Jk=0 V JeNftGTjj 



nf(k)>o/ 



«eyr,j *-0 \^JM,l^j, 



24. The method of claim 19 wherein the execution of step (h) is subject to 
the following constraint, 

25. The method of claim 19 wherein the execution of step (h) is subject to 
the following constraints. 



-A(r + l)=0 



A,(r)<i,, 

26. The method of claim 19 wherein the execution of step (h) is subject to 
10 the following constraint, 

yew, 

27. The method of claim 19 including the additional step (i) scheduling the 
movement of locomotives based on the result of step (h). 

28. An apparatus for optimizing the assignment of locomotives in a 
15 railroad network by minimizing the penalties associated with late train departures, 

wherein the penalty incurred, P^ , (r), if a train t at a location i does not have power 

over the predetermined time period r is known, said apparatus comprising: 

means for determining whether a train t is powered during the predetermined 
time interval r ; 
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means for assigning a first value to each train t if the train t does not have 
power and for assigning a second value to each train t if the train t does have power; 

means for multiplying the output from said means for assigning by the penalty 
incurred; 

means for summing output from said means for multiplying over all trains t 
moving between locations i and y, over all locations i, and over the predetermined 
time period r ; and 

means for minimizing the output from said means for summing. 

29. The apparatus of claim 28 wherein the means for minimizing is subject 
to the following constraint, 

30. The apparatus of claim 28 wherein the means for minimizing is subject 
to the following constraints, 

3 1 . The apparatus of claim 28 wherein the means for minimizing is subject 
to the following constraint, 

ZZauCt^)^^/. V.eF.re LIT,,. 

32. The apparatus of claim 28 wherein the means for minimizing is subject 
to the following constraints, 
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/eW, 

+iv:"(-)-er(-) 

- E a:;W-^-'(-K^r(--i)=TtEE^j:r(--^)' 



jtti 



t=0 ^yeW,(eT,^ 



33. The apparatus of claim 28 wherein the means for minimizing is subject 
to the following constraint, 

34. The apparatus of claim 28 wherein the means for minimizing is subject 
to the following constraints. 



E 



E«/.(^-6^y.-^y.J- E<.W 



-A(r + l)=0 



A(r)<X, 

35. The apparatus of claim 28 wherein the means for minimizing is subject 
to the following constraint, 

(?,,(r + lK^,,(r)>0 .V^eF,^e U T,,. 

36. The apparatus of claim 28 including means for scheduling the 
movement of the locomotives based on the output from -the means for minimizing, 

37. The apparatus of claim 28 including means for controlling movement 
of the locomotives in accord with the output from the means for minimizing. 

38. An apparatus for optimizing the assignment of locomotives in a 
railroad network by minimizing the cost of operating trains using a power class c, 
wherein the cost of using a locomotive of power class c as active power is known and 
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wherein the traveling distance of each train departing from location i is known, said 
apparatus comprising: 

means for determining the number of locomotives of power class c and 
signaling class s assigned to a train / initiated at location i over the predetermined time 
5 r ; 

means for calculating the product of the numb^ of locomotives, the cost of 
using a locomotive of power class c as active power in a train t and the distance train t 
will travel; 

means for summing the product over all power classes, over all signaling 
10 classes, over all trains traveling between locations i and j\ over all locations / on the 
network, and over the predetermined time t ; and 
means for minimizing the sum. 

39. The apparatus of claim 38 wherein the means for minimizing is subject 
to the following constraint, 

40. The apparatus of claim 38 whereui the means for minimising is subject 
to the following constraints, 

<5,,(r)<l VzeY.fe U T,, 

41. The apparatus of claim 38 wherein the means for minimizing is subject 
to the following constraint, 

42. The apparatus of claim 38 wherein the means for minimizing is subject 
25 to the following constraints. 
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JeNi JeN, 

nr(k)>of „ / \^ , X 



a 



nf{k)>of ^ 

43. The apparatus of claim 38 wherein the means for minimizing is subject 
to the following constraint, 

44. The apparatus of claim 38 wherein the means for minimizing is subject 
to the following constraints. 



A,(rKZ 



: Z<*{^-&jrOj- Z<(^) 



/6 JT, 



y«Ni 



~A(r + l) = 0 



A,(r)<A.. 

45. The apparatus of claun 38 wherein the means for minimizing is subject 
10 to the following constraint, 

46. The apparatus of claim 38 including means for scheduling the 
movement of locomotives based on the output from the means for minimizing. 

47. The apparatus of claim 38 including means for controlling the 
15 movement of locomotives based on the output from the means for minimizmg. 

48. A computer program for optimizing the assignment of locomotives in a 
railroad network by minimizing the cost of operating trains using a power class c, 
wherein the cost of using a locomotive of power class c as active power is known and 
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wherein the traveling distance of each train departing from location i is known, said 
computer program comprising the stq^s of: 

(a) determining the number of locomotives of power class c and signaling 
class s assigned to a train t initiated at location i over the predetermined time r ; 

(b) calculating the product of the result of step (a), the cost of using a 
locomotive of power class c as active power in a train t and the distance train t will 
travel; 

(c) summing the result of step (b) over all power classes; 

(d) sunaming the result of step (c) over all signaling classes; 

(e) suitmiing the result of step (d) over all trains traveling between 
locations z and j; 

(f) summing the result of step (e) over all locations i on the network; 

(g) svunming the result of step (f) over the predetermined time r ; and 

(h) minimizing the result of step (g). 
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